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Abstract 

We continue our study of the dynamics of the flavour sector of the Sakai-Sugimoto 
model in the presence of an external magnetic field, uncovering several features of 
the meson spectrum at high and low temperatures. We employ both analytical and 
numerical methods to study the coupled non-linear equations that result from the 
gravity dual. 
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1 Introduction 



Gauge/gravity dualities have become powerful tools for understanding diverse phenomena in 
strongly coupled gauge theories. In this paper we continue a line of investigation concerning 
the influence of an external magnetic field on the dynamics of the flavour sector of a gauge 
theory. The ultimate goal is to try to isolate features that are of interest for learning about 
the properties of nuclear matter, the idea being that these techniques provide an alternative 
(and/or complementary) route to that of Quantum Chromodynamics (QCD) for access to 
the physics of interest in various regimes. The model we study here has no pretensions 
to being close in microscopic detail to the known strong interaction physics (such as those 
captured by QCD), but nevertheless there is good reason to believe that there are phenomena 
that can be captured using it that are not sensitive to all of the detail. 

The model that we study here was, in essence, constructed by Witten in ref. |2] by taking 
the near-horizon limit of D4-branes compactified on a spatial circle to capture the physics of 
confinement/deconfinement transition in QCD. This model was later revisited in refs. [HE] 
incorporating the dynamics of probe D8/anti-D8 flavour branes, which led to a geometric 
understanding of spontaneous chiral symmetry breaking. It is now commonly known as the 
Sakai-Sugimoto model. 

In the Sakai-Sugimoto model the dynamics of the flavour sector has been analyzed at 
finite temperature in refs. [01 HH] and in the presence of an external electromagnetic field 
in refs. [TS| [T7] . For magnetic fields, it has previously been argued in e.g. ref. [H] that for 
(3 + 1) and (2 + l)-dimensional field theories there is a universal phenomenon at work that is 
sometimes referred to as magnetic catalysis of chiral symmetry breaking. These field theory 
treatments typically involve analyzing truncated Dyson-Schwinger equations in order to get 
access to the non-perturbative physics. 

This physics (that of a magnetic field catalyzing chiral symmetry breaking) can now be 
studied using the alternative technique of gauge/gravity duals (the magnetic field maps to 
a background B-field in the string theory and causes the flavour branes in the geometry to 
bend, generically) , and it has been observed [I8] to persist in all the Sakai-Sugimoto type 
models, where the flavour degrees of freedom are inserted as defects in the dual geometry for 
several (2 + 1)- and (3 + 1) -dimensional gauge theories. A similar effect has been observed 
to exist in other models in refs. [32l [36] . In those works, some of the structure of the meson 
spectrum was also studied, revealing phenomena such as Zeeman splitting, level mixing, etc. 
Here we would like to study the meson spectrum in the presence of an external magnetic 
field in the Sakai-Sugimoto model supplementing the study of the phase diagram obtained 
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in our earlier workllSI. 



Previous studies of the meson spectrum in this particular model have been carried out 
in e.g., ref. PS] for both low spin and high spin mesons. We uncover the role of the 
magnetic field in these spectra. For the low spin mesons we restrict ourselves to the high 
temperature phase only and therefore study the quasinormal modes of the scalar and a 
subset of the vector mesons. This comes with a caveat. The precise identification between 
the quasinormal frequency associated with the supergravity fluctuations in this model and 
the meson spectrum in the dual field theory is rather subtle and poorly understooci^. In 
the Sakai-Sugimoto type models the mesons transform under an t/(iVj)diag, whereas the 
quasinormal modes transform under either an U{Nf) l or an U{Nf)ji. This fact alone readily 
makes it difficult to precisely define what gauge theory quantities we would be studying 
by extracting the quasinormal modes. Nevertheless the quasinormal modes are naturally 
associated to the embeddings falling into the black hole of the gravity dual and therefore 
we study those in their own rights. 

Moreover, since the embedding function of the probe D8-brane is trivial {tq{u) = const, 
where tq{u) denotes the profile function) in the high temperature phase, the scalar and the 
vector fiuctuations remain decoupled even in the presence of a background anti-symmetric 
field such as a constant magnetic field. This is to be contrasted with the D3-D7 model 
in a background magnetic field [37] where the profile function of the probe D7-brane is 
non-trivial and therefore the scalar fiuctuations couple with the vector fiuctuations. Thus 
studying quasinormal modes in presence of a background anti-symmetric field in general is 
more difficult in the later model. 

The study of high spin mesons is another avenue that we pursue. We find that the 
presence of the magnetic field enhances the stability of such mesons by increasing their 
angular momentum. This is turn increases the dissociation temperature at which they fall 
apart into their constituents. We also uncover a simple realization of the Zeeman effect in 
the presence of a background field. It can be readily seen that such effects exist in all the 
Sakai-Sugimoto type models with probe branes of diverse dimensions. 

This article is divided into different sections. In section 2, we review the Sakai-Sugimoto 
model and its phase diagram in presence of a magnetic field. We then move on to section 
3 to study the quasinormal modes for scalars and the decoupled sector of transverse and 
longitudinal vector mesons in this model. Section 4 is devoted to the study of high spin 
mesons, where we first obtain some analytical results and then proceed to the numerical 
analysis. We have relegated the details of some of the calculations to four appendices. 



^For a detailed analysis regarding the issue of precise identification, see ref. [3T] 
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Appendix A contains the details of the Lagrangian for the quadratic fluctuations and the 
equations of motion for the gauge fields as well as a detailed discussion about when they 
can be decoupled. In appendix B we study the longitudinal and transverse modes in the 
hydrodynamic limit in vanishing magnetic field and obtain the well-known hydrodynamic 
dispersion relation for the longitudinal mode. Appendix C contains the general variable 
changes we use to recast the equations of motion for the fluctuation modes into an effective 
Schrodinger equation. In appendix D we present a model calculation of high spin mesons 
in a background magnetic field by analyzing spinning strings in Rindler space coupled to 
a magnetic field. This serves as a toy model that captures some of the features of meson 
dissociation in a background field. 

2 Review of the Model 

The Sakai-Sugimoto model is constructed from the near horizon limit of D4-branes with 
one compact spatial direction along which an anti-periodic boundary condition is imposed for 
the adjoint fermions. This makes the adjoint fermions massive and breaks supersymmetry. 
Flavours are introduced by Nf D8 and Z^S-branes intersecting the D4-brane at a (3 + 1)- 
dimensional defect. Thus there is a global U{Nf)L x U{Nf)fi flavour symmetry as viewed 
from the worldvolume of the Z}4-brane. We work in the probe limit, namely Nf <^ where 
the flavour branes do not backreact on the background. 

At low temperature the near horizon geometry of the D4-brane is given by 

where are the flat 3-directions, is the Euclidean time coordinate, r is the spatial 
compact circle, dfll is the metric on the round and u is the radial direction, ig is the 
string length, gs is the string coupling; V4 and 64 are the volume and volume form of 
respectively. Also, cj) is the dilaton and F(4) is the RR four-form field strength. The radius 
of the compact spatial circle is given by 27ri?4 = (47ri?^)/(3M^^), which is obtained from 
the condition that the compact direction r shrinks away smoothly ai u — ukk- 

The five-dimensional gauge coupling is given by g1 = {27Tygsis, the four dimensional 
gauge coupling can be obtained by dimensional reduction yielding gf = g'^/27TR and the 
five dimensional 't Hooft coupling is defined to be A = {g'^Nc)/4:'K. Gravity calculations 
are reliable when the spacetime curvature is small compared to the string tension and also 
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when string loop effects are suppressed. Talcen together these two conditions tell us that 
the gravity approach is valid in strong 't Hooft coupling limit of the four- dimensional gauge 
theory (see e.g. ref. 0). 

The high temperature background is given by 

ds' = {^f [dxdx' + f{u)dtl + dr') + (-^ + uH^i^ , 

where T is the background temperature and all other parameters are as defined before. 

We let the probe brane-anti-brane pair stretch along the directions {t,x*,f24} and have 
a non trivial profile described by the function t{u). In the low temperature phase the 
background cigar topology in the {r, u}-submanifold forces the brane pair to join at some 
radial distance Uq breaking the global U^Nf)^ x U{Nf)ji down to a f/(A^/)diag- This is 
then identified with the spontaneous chiral symmetry breaking in the model. The high 
temperature phase is richer in content. The cigar topology now appears in the {t^;,^}- 
submanifold of the background and therefore the brane-anti-brane pair can end on the 
horizon separately. Depending on the background temperature or the asymptotic separation 
between the brane-anti-brane pair, there is a chiral symmetry restoring transition. 

The background geometry undergoes a confinement /deconfinement transition at a tem- 
perature Trf = 1 / {271 R4), where R4 is the radius of the compact r-direction. This is the 
temperature at which the gluons deconfine. When flavours are added, depending on the 
dimensionless parameter L/R4 there can be a range of temperatures where the gluons are 
deconfined and chiral symmetry in the flavour sector is spontaneously broken by virtue of 
the brane-anti-brane pair joining at the radial position Uq. If the quark separation L obeys 
the bound that L/R4 < 0.97, then the intermediate phase (which is deconfined but chiral 
symmetry broken) exists. If however L/Ri> 0.97, then deconfinement and chiral symmetry 
restoration happens together. Therefore there can exist a separation of scales between the 
gluon deconfinement and chiral symmetry restoring transitions. This has been discussed in 
refs. P[in]. 

We can also set i?4 — >■ cxo so that the spatial circle direction now becomes a flat extended 
direction; this particular limit has been studied in detail in refs. [6[[7l[8]. It turns out that in 
this particular limit the dual gauge theory is described by the Nambu-Jona-Lasinio model 
(when probing with D8-branes) or the Gross-Neveu model (when probing with D6-branes) 
with a non-local four-fermi interaction. In this case, the asymptotic separation between the 
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brane-anti-brane pair sets the coupling for the four-fermi interaction. Here the so called 
intermediate temperature phase always exists. 

T 
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Figure 1: The phase diagram vaT — B plane. The vertical axis is expressed in units of 1/L and 
the horizontal axis is expressed in units of R. 

In this model, the effect of a background magnetic field on chiral symmetry breaking was 
analyzed in refs. [171 [IE]- It was explicitly shown that a background magnetic field enhances 
chiral symmetry breaking and increases the chiral symmetry restoring temperature. The 
resulting phase diagram is shown in figure [H This magnetic field couples to the flavour 
sector only, therefore it does not affect the (gluon) confinement /deconfinement transition 
temperature. Thus the background magnetic field increases the separation of gluon decon- 
finement and chiral symmetry restoring temperature scales. However, as we can observe 
from figure [1] this separation of scales does not increase indefinitely with increasing mag- 
netic field, rather it tends to saturate a maximum value. This particular feature persists in 
several other Sakai-Sugimoto type mo dels p^. 

3 Mesons with Small Spin 

Mesons with low spin in the dual gauge theory correspond to the small fluctuations of 
the classical profile of the probe brane. The fluctuation in the geometric shape (which is 
denoted by the function t{u)) of the probe branes corresponds to the scalar (or pseudo- 
scalar) meson mode and the fluctuation of the Maxwell field on the worldvolume of the 
i58/D8-brane corresponds to the vector (or axial vector) mesons. 
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3.1 Scalar Fluctuation 

Let us consider a fluctuation of tlie probe brane embedding function t{u) = to{u)+x, where 
tq{u) is the classical profile. Before proceeding further let us recall that the function tq{u) is 
non-trivial only for the U-shaped configuration (for the straight branes, this function is just 
a constant). Here we will consider the high temperature phase only for which tq{u) = const, 
and the scalar fluctuation is always decoupled from the vector fluctuations. 

To recast the equation of motion for the fluctuation mode into an equivalent Schrodinger 
equation we define, fx{u) = a{u)g{u), where a{u) is defined in eqn. (IC-2|) (the details of the 
coordinate changes are given in appendix C). The equivalent Schrodinger equation is given 
by (to avoid notational clutter we choose R = 1 and Ut = 1) 



Here u is the "tortoise" coordinate defined in appendix C. For simplicity, we continue to 
express the Schrodinger potential in the original u-variable. 

We consider only in-falling modes at the horizon; such a solution can be written as g{u) = 
exp{—iujt)tp{u) . Now multiplying the first equation in ([3]) by iIj*{u), integrating by parts 
from the horizon to the boundary and using the equation of motion we get 



Now we can understand qualitative features of imaginary part of quasinormal frequency 
by studying this effective potential. We resort to numerics to study the potential K(w). 
This potential vanishes at m = 1 because of the presence of the horizon and diverges as 
u —>■ oo reflecting infinite gravitational potential of the background geometry. 

The important feature that we find from figure [2] is that the Schrodinger potential is 
always positive, therefore from eqn. (jl]) we can conclude that the imaginary part of the 
quasinormal frequency is strictly negative. Moreover, the potential does not develop a 
negative energy well, therefore there exists no negative energy bound states in this phase. 
This is consistent with the fact that the probe brane configuration is stable and tachyon free 
and chiral symmetry restoration is accompanied by meson melting transition. Tachyonic 
modes in the scalar field fluctuation are likely to appear in the overheated low temperature 
phase which is bypassed by the transition. Also we see that the potential vanishes at the 
horizon and therefore the spectrum is continuous, which is consistent with the melting 
transition. 




and Vs{u) 



{u^ - 1) (554 {7u^ + 5) + 2B^u^ (71^3 + 7) + 16u^ {5u^ + 1)) 

16u5 (52 + u^f 



(3) 
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Figure 2: The Schrodinger potential for different values of external magnetic field. The dashed 
vertical black line represents the position of the horizon, the top-most curve represents the po- 
tential at zero external field. As the field increases we find the potential grows more slowly with 
increasing u. 



Let us now represent the quasinormal frequency as = cur + iujj. It can be noted from 
the effective Schrodinger equation that there is a Z2 symmetry in cur (namely, if g{u) is 
a solution with a frequency 00 then g*{u) is also a solution with frequency — a;*[ll]). Now 
we can observe that the Schrodinger potential increases monotonically as we move away 
from the horizon. This implies that Re[u;^] > and hence Abs[u;R] > Abs[ct;i]. By explicit 



numerical calculation we will find that this is indeed the case (we point to figure 3(a) and 



3(b)). 



To extract more of the physics, we now compute the quasinormal frequency corresponding 
to this scalar fluctuation. We find it extremely convenient to perform a change of coordinates 
as implemented in ref. [30]. Let us introduce a new variable p = u/Ut and v = t + a{p), 
where a{p) is determined from the condition that 

dv = dt+ „]. . dp . (5) 

pv(p) 

This coordinate system {v, p} makes the numerics conveniently stablj^ to find the quasi- 
norma, modes. The equation of .notiorfi for the flnctuation now takes the form 

d, {p''/\l + B'p-'f'fr^) - luj [d,p'l\l + B'p-'f'f^ + 2p'l\l + B^p-'f/'f^] 

V/'(l + 5V')i/Vx = 0. (6) 



^We refer the reader to ref. [30] for more details. 

^In presence of an external field this equation does not take the general form described in ref. [30]. This 
is due to the breaking of the 50(3, 1) S'0(1, 1) x SO{i) by the presence of the magnetic field as a result 
of which the Laplacian in the {x^, a;"^}-direction gets squashed. 
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It is now straightforward to check that this equation admits the existence of quasinormal 
modes. To find how it depends on the background magnetic field we use a shooting technique 
to solve the fiuctuation equation and pick the appropriate value of uj for which the equation 
admits a normalizable solution (this is achieved by imposing the condition that — as 
p oo). As boundary conditions we impose fy^{l + e) = 1 and /^(l + e) equal to a value 
obtained by requiring the equation of motion to be regular near p = 1, the event horizon. 
Here the parameter e is a vanishingly small number in our numerical scheme. 
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Figure 3: The real (wr = Re[a;]) and imaginary {uji = Im[a;]) part of the quasinormal modes for 
the scalar fluctuation as a function of the external field. The quasinormal frequency is measured in 
units of the background temperature T. The small wiggles are due to numerical errors; here we have 
plotted the polynomial which fits the data. As expected before, we confirm that Abs[wR] > Abs[a;i], 
however they are roughly of the same order. 



The result of this pursuit is summarised in figure 3(a) and 3(b) Since we are in the high 
temperature phase the quasinormal mode is comprised of both real and imaginary parts. 
It would be expected that the absolute value of the real part corresponds to the mass of 
the meson before it melts and the absolute value of the imaginary part corresponds to its 



inverse lifetime analogous to ref. [H] for the D3-D7 system. Figure 3(a) then tells us that 



the magnetic field increases the scalar meson mass and figure 3(b) hints that the inverse 
lifetime gets bigger as the field is dialed up. However the precise identification of the meson 
mass and inverse lifetime with the quasinormal modes of the scalar fiuctuation is a subtle 
issue in this particular model (see ref. [31]). 



3.2 Vector Fluctuation 



The vector meson spectrum in general is determined by a set of coupled equations of motion. 
It is however possible to completely decouple the transverse and the longitudinal modes if 
we restrict ourselves to the oscillations parallel to the magnetic field. If we restrict ourselves 
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to the oscillation in one of the perpendicular directions (e.g. along the x^-direction) to the 
applied magnetic field, then it is also possible to decouple one of the transverse modes (in 
this case the 742-vector mode). However, the remaining transverse mode and the longitudinal 
mode remain coupled. The details are given in appendix A. Here we study the decoupled 
sectors only. 



3.2.1 The Transverse Mode 



We begin analyzing the transverse vector meson spctra. To extract qualitative features of 
the spectra we study the Schrodinger potential. The equation of motion for the A2 vector 
modes is given by 



du 



- det (Ei;))5^^5"M',(«) + e- det (eIJ)) co^S^'S^'Mu) = , (7) 



where we chose an ansatz of the form A2 = A2{u)exp{—iujt) . It is straightforward to turn 
this equation in the form of a Schrodinger equation with an effective potential given by 

{u^ - 1) (5 (7^3 + 5) 5^ - 2 {u^ - 43^3) 5^ + 8 {u^ + 2u^)) 



V.u) 



16^5 (m3 + B^f 



^3/2 



with du = — du 

- 1 



where u is the "tortoise" coordinate. Clearly the horizon is located at m — >■ 00 and the 
boundary is located at m — 0. 

The corresponding potential is shown in figure |H The main qualitative feature that 
emerges from this is the fact that there is no instability due to existence of bound states 
of positive imaginary frequency modes since the potential remains positive definite. We 
observe, however, that a non-vanishing magnetic field can create a dimple in the otherwise 
smooth and rather featureless potential and therefore create local (positive) potential well. 
Therefore it is possible to have long-lived bound states with Abs[a;R] ^ Abs[co'i]. In this 
case, Abs[a;R,] is set by the depth of the well, which in turn is set by the background field. 
If we increase the magnetic field even higher then the dimple gets bigger and bigger and 
ultimately swallows the well. At even higher magnetic field the potential again takes a 
shape similar to the vanishing external field case. Therefore the imaginary part of the 
quasinormal frequency may have an interesting behaviour with increasing magnetic field. 
We may expect that Abs[c(Ji] starts decreasing with increasing magnetic field, but after a 
critical value it starts increasing again. Next we turn to the numerics. 
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Figure 4: The Schrodinger potential for the vector meson for different values of external magnetic 
field. Here u denotes the "tortoise" coordinate mentioned in appendix B. As the magnetic field 
increases we observe that a small well develops in the effective Schrodinger potential. The depth 
of this positive well is set by the background field. 



The equation of motion (for the vector mode A2 with non-zero spatial momentum) in this 
case is given by 



da 



0, ZG{1,3} 



(9) 



where S""^ are the analogue of S"'^ in eqn. flA-4p when the coordinate change described in 
eqn. ([5]) has been performed. 

The background field breaks the full SO (3, 1) Lorentz symmetry and there is an unbroken 
SO (2) symmetry corresponding to the rotation in the plane perpendicular to the magnetic 
field. Therefore the transverse vector meson spectra with a non-zero momentum will have 
two distinct branches. Mesons can have momentum along the direction of the background 
field (corresponding to i = 1) or in the perpendicular plane (corresponding to i = 3). 

Without any loss of generality we can consider studying the spectrum of the A2 vector 
meson. We can give this meson a momentum along the x^-direction or along the x'^-direction 
(but no momentum along x^-direction since this is the transverse mode). The spectrum 
would be entirely equivalent to the spectrum of the vector modes having momentum 
along the x^-direction or along the x^-direction. 

Now we follow the same numerical approach. We impose ^2(1 + e) = 1 and fix ^2(1 + e) 
from the equation of motion near the horizon. With these boundary conditions we look for 
normalizable solutions for ^2, which gives the quasinormal modes. 



In figure 5(a) and 5(b) we have shown the dependence of the real and imaginary part 
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Figure 5: The real and imaginary part of the quasinormal modes for the vector fluctuation as 
a function of the external field. The frequency is expressed in units of background temperature 
T and the magnetic field is expressed in units of R. We find a local minima in Abs[u;i] in the 
parameter range we have explored. This is qualitatively similar to the behaviour anticipated from 
analyzing the effective Schrodinger potential. We also observe that in the most of the parameter 
space Abs[LL'R] » Abs[a;i]. 



of quasinormal frequency for the particular choice of A; = 0. We observe that the inverse 
lifetime of the vector meson gets shorter as the external field is dialed up and therefore 
enhances the stability of the meson. This behaviour is opposite to what we observed for the 
scalar meson fluctuation. 



Next we obtain the dispersion relation when k ^ 0. In figure 6(a) we demonstrate this 
dispersion relation for a wide range of values for k and in figure 6(b) we have shown a 
magnified version of figure 6(a) to more clearly see the role an external magnetic field plays 
in this case. These two figures correspond to the dispersion relation when the momentum 
is parallel (i.e. along x^-direction) to the magnetic field. We observe that for high enough 
spatial momentum the dispersion relation curve approximates the zero field dispersion curve 



quite well. We have shown in figure 6(c) the dispersion relation when the momentum is 
perpendicular (i.e. along x^-direction) to the background field. The qualitative behaviour 
between these two dispersion curves are clearly different; however for large enough spatial 
momenta these curves tend to become insensitive to the background field. This is because 
the heavier the meson becomes the less sensitive it is to the background field. 
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Figure 6: The dispersion relation between the real and imaginary part of the quasinormal fre- 
quency for two different values of magnetic field. The red (solid) curves correspond to the disper- 
sion relation for non-zero magnetic field. The black (dashed) curve corresponds to the dispersion 
relation when there is no magnetic field. The black dots represent the corresponding values of 
quasinormal frequency for vanishing spatial momentum. Figure 6(a) and |6(b)"] correspond to vec- 
tor mesons having momenta along x^-direction and figure 6(c) corresponds to vector meson with 
momentum along x^-direction. 



3.2.2 The Longitudinal Modes 



Now we study the longitudinal mode. With the change in variables in equation ([5]), we write 
down the equations of motion for gauge fields along {t>, x^} 



e-M-det (^eS!) {S^'^'Spp - S'^'S^") A'^ 



-det (^1;^)A; [^™(fc^ + LuA,) + z^^M'J + Agj'BA', = 



e-'^J-det {SPPA[ - iSf"'{kA, + luAi)) 



-e-M-det (^E^^)s'^iu [5™(M, + cuAi) + tS'^'A',] - Ag^^BA', = , (10) 
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along with the constraint 



e-'^W -det (eIJ) {SppS'''' - S'^'S'"') tuA'^ + t4g;^B{ujA^ + kA,) 



-ge-M-det (^1?)^'' [iS^pA', + SP'^ikA, + luAi)] = , 



where k is the momentum along x^-direction. We define the longitudinal mode to be given 
by 

uAi + kA^ = E . (12) 

Now using the definition of the longitudinal mode and the constraint equation we can solve 
for A'^ and A'^ in terms of E and E' and then substitute back in the equation of motion for 
the gauge fields to obtain 



E"(p) + Zi(p)E'(p) + Z2{p)E{p) = 



(13) 

where Zi and Z2 are known function^ of p. 

Equation (fT3l) now can be numerically solved by using the boundary conditions that 
E(l + e) = 1 and E'(l + e) equal to a value obtained from the equation of motion itself. The 
condition of normalizability then fixes the quasinormal modes. 
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Figure 7: The dispersion relation of longitudinal oscillation in presence of a background magnetic 
field. The quasinormal frequency is measured in units of the background temperature. In figure 



7(a) the black solid curve represents the dispersion relation when the magnetic field is absent, 
the blue (circular) dots correspond to i? = 5 x 10~^ and the purple (square) dots correspond to 
B = 10~^. B has been expressed in units of R. In figure [7(b) | we have shown the dependence 
of the lowest quasinormal frequency corresponding to the hydrodynamic mode with the magnetic 
field for a fixed momentum k = 0.01. The imaginary part of the quasinormal frequency does not 
depend on the magnetic field in the hydrodynamic limit. 



"^The explicit expressions of Zi and Z2 are not very illuminating, thus we do not provide them here. 
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The numerical result is showij^ in figure 7(a) The qualitative behaviour of the dispersion 
curve is similar to the transverse modes. It is evident that increasing magnetic field brings 
quasinormal frequency closer to the real axis in the complex plane. However, for a fixed 
magnetic field, the dispersion curves tend to approach the zero magnetic field dispersion 
curve as the momentum increases. 



4 Mesons with Large Spin 

Mesons with large spin can be described by strings having both its end points on the flavour 
brane and rotating on a two plane along the directions where the gauge theory lives. For 
large enough angular momentum, the length of the string is much larger compared to is and 
we can use classical Nambu-Goto action to describe the dynamics of the rotating string. 
Mesons with large spin have previously been analyzed in, e.g., refs. [25l [20] . Here we will 
follow the same framework, now including the effect of an external magnetic fleld. 

In absence of magnetic fleld due to rotational symmetry, the spectrum does not depend on 
the direction of the angular momentum J. An external magnetic fleld breaks this symmetry 
and introduces a splitting of the energy levels (Zeeman effect). Here we will focus on the two 
cases, namely when the magnetic fleld is parallel or anti-parallel to the angular momentum. 
In this case there should be a Zeeman splitting in energy for mesons with non-zero magnetic 
moment. 

This idea has been implemented in ref. [H] in a different modeH We adopt similar system 
and consider two flavour brane-anti-brane pair very close to each other so that in the bulk 
they join leaving two diagonal f/(l)s (call this U{1)a and U{\)b corresponding to the two 
pair of flavour branes). We consider now turning on flelds having equal opposite charge 
under these two t/(l)s. Now AA and BB strings will have equal and opposite charges but 
AB and BA strings will have equal charges. In the flrst case, the total orbital magnetic 
moment of the meson will vanish whereas in the second case the meson will have a non-zero 
magnetic moment. The magnetic fleld will induce a quadratic correction to the meson energy 
for the former and a linear Zeeman splitting for the later. Here we will analyze the later 
case only. It is useful to note that this conflguration is symmetric under reflection around 
the midpoint of the string (at p = 0). We will refer to this as the symmetric conflguration. 

^Note that here we numerically explore only for positive values of B. The equation of motion for the 
longitudinal modes have a linear term in B, and therefore the quasinormal modes may depend on the sign 
of B. However, we found numerically that the qualitative features of the quasinormal modes are similar 
also for negative values of B. 

^For more recent similar work see, e.g., refs. [421143] , 
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To proceed we write the relevant part of the background metric^ in the following form 
(rewriting the {x'^,x^}-two plane in polar coordinate {p, 0}) 



1i \3/2 

'r 



) [f\u)dtl + dp" + p'^dcP^) + [-\ 



du' 



(14) 



We recover the zero temperature background by setting ut = ^- We will work with the 
Nambu-Goto action for the string with the following ansatz for the string profile 



t = T , p = p{a) , u = u{a) , = UJT . 



(15) 



From now on the parameters {t, cr} will always refer to the worldsheet coordinates for the 
string. We also assume that u is positive (i.e. only clockwise rotation). 

In the presence of an external constant magnetic field the Nambu-Goto action is accom- 
panied by a boundary term coupled to the magnetic potential = Bp'^/2. This potential 
gives rise to a magnetic field 5(2) = Bpdp A dcf) = Bdx^ A dx^. Taking this boundary 
term into account the action for the string is given by 



S 



AS 



27Ta' 



drda 



B 



2Txa' 




1/2 



B 



(16) 



where a represents the right and the left boundaries and the relative positive sign between 
the two boundary contributions coming from the magnetic potential is due to the symmetric 
configuratioij^. For convenience we introduce the rescaled variables 



p = pcu , u 



uuo 



R = Ruj 



:i7) 



The expression for the energy and the angular momentum for the spinning string in the 



^Here we explicitly write down the metric corresponding to the high temperature phase only. The 
background metric corresponding to the low temperature phase can likewise be written in an analogous 
form. 

^For the non-symmetric configuration the there will be a relative negative sign between these two terms. 
We refer to ref. [41J for further details on this construction. 
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rescaled variables can simply be obtained to give 
dL ^ 1/1 



E 



ijj- L 



Ilia' 



2Txa' \ ujy^R^I^ ) 



^5/2^3/2 



~/2 I ~/2 
P +U 



da- 



p2 




(4 


.*4 







(18) 



We observe that the external field can enhance angular momentum of the meson if it is 
aligned in parallel (positive values of -B); and can reduce the angular momentum if it is 
anti-parallel (negative values of B). 

To proceed we need to solve for the string embeddings. The equation of motion obtained 
from the Nambu-Goto action has to be supplemented by the boundary condition 

dL 



TT 



,5X 



M 



. 



(19) 



9S 



Since 5u\qj. = and Splgj. is arbitrary, we have to impose the Neumann boundary condition 
(in absence of any external field) {dL/dp')\Q^ = 0. This condition gives us the following 
constraint 



TT 



11 

pies 



P 



{u/Rwyii-p 

P + m'2 {ujR/u) 



1/2 



. 



(20) 



9S 



So we need to impose the boundary condition p'\qy, = O5 which means that the string is 
ending perpendicularly on the fiavour brane-anti-brane pair. We will come back to the 
details of our numerical scheme in a later section. 

The presence of an external magnetic field changes this boundary condition. It can be 
easily seen that the boundary condition to be satisfied is now given by 



vr 



.11 

p\dT. 



-Bp) = 



(21) 



where ^rlj^^ is given in equation (!20|) . Therefore we need to solve the equation of motion 
obtained from the action in eqn. fll6p subject to the boundary condition given in eqn. fl20p 
or fl2T|) in absence or presence of a magnetic field respectively. 

We will begin by obtaining some analytical results in the large angular frequency limit in 
the zero temperature phase. We will choose a = p gauge and send uj —>■ 00; from the reality 
condition of the Nambu-Goto action in eqn. (fT6l) this limit sets a bound \p\ < l/u ^ 0; 
and therefore equivalently we study the short string limit. 
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4.1 Analytical Results 



In this section we consider the zero temperature background and spinning strings with large 
angular frequency only. 



4.1.1 Vanishing Magnetic Field 

In this section we closely follow the approach adopted in ref. [20]. In the large angular 
frequency limit the relevant string profile is well approximated by the local geometry very 
close to the point from where the string hangs. In this limit the string hangs from the point 
where the brane-anti-brane pair join (the radial position denoted as Uq). We therefore pick 
the following ansatz to approximate the string profile 

u(p) =u;Uo + ^ . (22) 

00 

It is easy to check that this ansatz is consistent with the equation of motion for the string 
in — i> oo limit. To obtain the leading order behaviour we have to solve for the function 
/(p) with the boundary conditions that /(O) = = /'(0)§. Substituting the ansatz given in 
eqn. fl22l) in the equation of motion for the string and keeping only the leading order terms 
we obtain a differential equation for the function /(p) which can be integrated to obtain the 
following analytic form 

/(P) = |^(p' + arcsin^(p)) • (23) 
Plugging back f{p) from eqn. ( l23i) to eqn. ( fTSi) and simplifying we get the following relations 

From this we can see that uj ^ oo limit is equivalent to J ^ A limit. Also we find that the 
mesons with large angular frequency follow a Regge trajectory with an effective tensioijl° 

Teff = EV(27rJ) = Af/o'/V(2i?'/'). (25) 



4.1.2 Non-vanishing Magnetic Field 

In presence of a magnetic field the leading order solution for u{p) remains the same as in eqn. 
(l23ll : however the boundary condition needs to be modified as given by eqn. (1211) . Using 

^At u — loUq we require u'{p) = 0, which sets the boundary condition for the function f{p). 
^"The radial position Uq is related to the asymptotic separation L between the flavours via J7q^^ = 
E?/"^ / (4L)i3(9/16, 1/2), where -B(a, 6) is the Beta function. This length L sets the coupling strength of the 
non-local four Fermi interaction for the dual NJL model. 
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the expression for vr^l^^ from eqn. (!20l) (in the gauge cr = p) in the modified boundary- 
condition given in eqn. fl2T|) we can determine the value of p at the boundary (meaning at 
u = Uq where the string ends on the flavour brane). In the limit of small magnetic field this 
critical value of p (denoted as pc) can be obtained to be 

Pc = l-^ TT • (26) 



With this pci there will be two different string profiles which we have shown in figure 8(a 



and 8(b) respectively. 




u(o-) 



u(o-) 




P(o-) 




P(o-) 



(a) 



(b) 



Figure 8: The two possible profiles are presented; Ug is the radial turnaround position for the 
string. These are obtained by solving the equations of motion numerically, although we come back 
to the numerical results in a later section. The profile in the left hand side has more energy than 
the one in the right hand side. 



These two different profiles result in two different energies (let us denote these energies by 
E±). The one with the lower energy (i.e., with ES) is extended from p = to pc and the 
one with the higher energy (i.e., with Ej^^ is extended from p = to p = 1 and then folds 
back to p = Pc- Using the formulae for meson energy E and angular momentum J given 
in eqn. (ITSl) we get the following linear order corrections (Zeeman splitting) introduced by 
the magnetic fielci^ 

E. = y^^vxj ^ 



^+ - ^^'^^^ + 7^Va-^- ^^^^ 



"'^^The parameter Uq is also related to the constituent quark mass. For a study of how the constituent 
quark mass behaves (for a fixed value of L) with the magnetic field see ref. [18] . 
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It is easy to check that this sphtting is consistent with the general formula for energy 
(obtained in ref. [lO]) of a Nambu-Goto string with an effective tension Tcs moving in a 
background magnetic field. The energy spectrum (for large J) is given by the following 
formula 



Recalling the expression for Tes from eqn. (l25l) and using the general result above the linear 
order correction to the energy of the string can be obtained which matches exactly with 
what is obtained in eqn. fl27|) . 

4.2 Numerical analysis 

In order to extract the relevant physics we need to solve for the string embedding. To that 
end we choose the parametrization a = p + u since it provides better numerical stability as 
discussed in e.g., ref. [H]. 

To solve the equation of motion first we fix the maximum radial position the string can 
attain (this is the position where the brane-anti-brane pair join, which is denoted by Uq). 
This in turn fixes the asymptotic separation between the brane-anti-brane pair. Fixing 
Uq is equivalent to fixing the constituent quark mass of the dual gauge theory and fixing 
L, the asymptotic separation of the brane-anti-brane pair, is equivalent to keeping the four 
fermi coupling of the dual NJL model. Now for a given Uq we look for a radial turnaround 
position for the string (call it Ug) such that shooting from Ug with the IR boundary conditions 
u{(7)\fj^ = Us and u'{a)\jj^ = would satisfy the UV boundary condition M'(cr)|^^^ = 1 (up 
to numerical accuracy). The existence of the turnaround position is guaranteed by the 
symmetric configuration that we consider here. Below we show two such representative 
profiles for u = 0.5 and u = 3 respectively (corresponding to the blue and the red curve). 

Such profiles were discussed in ref. [20] in a different but closely related model. In brief, 
we find that for a given value of u there are classes of solutions distinguished by the number 
of times they cross zero (referred to as "nodes" in ref. [20]) along the horizontal axis. Here 
we will constrain ourselves to the n = 1 case only. 

We will now explore the behaviour of high spin mesons in presence of a background 
magnetic field. The key feature that we will observe for both zero and finite temperature 
phases is the existence of a shift in physical quantities such as energy and angular momentum 
in presence of a non-zero magnetic field. We begin with our numerical results for the zero 
temperature phase. 
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U((T) 






Figure 9: Two representative spinning string profiles. We have set R = 1. 



4.2.1 Zero Temperature Phase 




0.5 1 1.5 2 2.5 3 



(c) 

Figure 10: The effect of an external magnetic field in the zero temperature phase. The black 
dashed lines represent i? = 0, blue dashed lines represent B = —0.1 and the red dashed lines 
represent B = 0.1; E and J have been evaluated in units of l/(27ra'). B have been evaluated in 
units of R. 
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In figure 10(b) we have shown the dependence of the meson energy E a.s a function of 



the angular frequency u, in figure 10(c) we have plotted the dependence of meson angular 



momentum again as a function of u and in figure 10(a) we have shown the dependence of 



energy with angular momentum of the meson. For zero external magnetic field we recover 
the behaviour discussed in ref. [25] . We find that the external magnetic field shifts the 
energy and the angular momentum of the meson for generic values of the parameters of 
the system; however with increasing angular frequency these observables tend to become 
insensitive to the external field. In view of the approximate analytic results obtained in eqn. 
(1271) it is straightforward to observe that an external field contributes only at order B/u 
and hence has vanishingly small effects for large enough angular frequency. 

The existence of two energy branches for the same given value of J seems to persist 
(as in vanishing external field) even in the presence of a magnetic field, although as the 
magnetic field is increased we observe a tendency that this multi-valuedness in energy starts 



disappearing. Therefore any possible unstable upper branch in figure 10(a) can get promoted 
to a stable one by having sufficiently high angular momentum. 

Moreover we observe that there exists a maximum angular momentum beyond which the 
spinning meson dissociates. This dissociation is mediated by the acceleration of the spinning 
string. However the presence of magnetic field can again stabilize the mesons by raising the 



maxima in figure 10(c) Physically this simply corresponds to the fact that the magnetic 



field enhances the angular momentum of the meson and thus makes it stable. 




B=0 



10 
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Figure 11: Dependence of maximum spin with magnetic field for varying constituent quark mass. 
Jcrit and rriq have been calculated in units of l/{2ira'). We observe that even for a certain negative 
value of B, i.e. when the magnetic field is aligned anti-parallel to its angular momentum, the 
meson can be more stable than its zero field counterpart. 



To demonstrate this effect we observe the plot in figure [TTl Clearly the B = —0.1 branch 
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is higher than the B = branch; an even higher value of B would inhibit any dissociation 
at all. Also the heavier the constituent quark mass is the more Jcrit becomes making it less 
likely for the heavy mesons to dissociate. 

E 



"""-*--«., 8.2 




^^•-^ 8.1 








7.9 
7.8 




7.7 





B 



-1.0 -0.5 0.0 0.5 1.0 



Figure 12: Linear correction to energy for the equal charge case. From top to bottom the values 
of angular momenta are J = 0.937, 0.907, 0.862 respectively. The dashed straight lines are the 
best fit lines. E and J are both presented in units of l/(27ra'). We observe that all the three 
best fit lines have same (up to numerical accuracy) slope, which we expect from the very nature 
of Zeeman splitting. 



A key feature of the high spin meson spectrum in presence of an external magnetic field 
is the linear Zeeman splitting. In figure [12] we have shown such linear correction in energy 
for fixed values of the meson angular momentum. 

4.2.2 Finite Temperature Phase 

We now perform a similar numerical study of high spin mesons in this intermediate tem- 
perature phase. 



In figure 13(b) we have shown the dependence of meson energy as a function of its angular 



frequency oo; in figure 13(c) we have plotted the angular momentum of the meson as a 



function of uj and in figure 13(a) we have shown the dispersion relation between the meson 
energy and its angular momentum. The qualitative features of finite temperature physics 
are similar to that of the zero temperature physics. There is however one distinction as 
compared to the energy spectrum at zero temperature phase. The presence of magnetic 
field at finite temperature does not promote the possible unstable upper branch in figure 



13(a) to a stable one within the range of values for the magnetic field that we have explored 



using our numerical approach (we expect however that for sufficiently high values of B this 



23 





Figure 13: The shifts in presence of an external magnetic field. E and J have been computed in 
units of l/(27ra'). B has been evaluated in units of R. The black curve corresponds to vanishing 
magnetic field, the blue curve corresponds to setting B = —1 and the red curve corresponds to 
setting B = 1, where B has been expressed in units of R. 



unstable mode will be promoted to be a stable one). A possible thermally enhanced decay 
channel therefore remains open at finite temperature for a rather high value of the external 
field. The other familiar role of the magnetic field, we again find, is to introduce shifts in 
the physical quantities such as the energy and angular momentum of the meson. 

The meson dissociation due to spin should now be enhanced due to thermal fluctuation. 
We have previously observed that the magnetic field adds extra angular momentum to the 
system and stabilizes the mesons. We study in figure 14(a) the effect of two competing 
parameters on the critical angular momentum of the meson. Finally we observe the familiar 



Zeeman splitting in figure 14(b) 
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•^crit 




(a) (b) 



Figure 14: Figure 14(a) shows the dependence of maximum spin calculated in units of l/(27ra') 
with the applied magnetic field. We can observe that there exists a certain range of values for 
Uq, which in turn fixes the value of the constituent quark mass, where the angular momentum is 
lowered by the magnetic field. Figure 14(b) shows the linear correction to meson energy for the 
equal charge case at finite temperature. From top to bottom the values of angular momenta are 
J = 0.630, 0.560 and 0.490 evaluated in units of l/(27ra') respectively. The dashed lines are the 
best fit lines. 



5 Conclusions 



We have uncovered a sub-sector of the meson spectrum in the presence of an external 
magnetic field, complementing the phase structure that we obtained in ref. [18]. Mesons 
with small spin that can be obtained by studying the quadratic fluctuations of the probe 
brane classical configuration are harder to study in absolute generality. In general the task is 
to solve a set of coupled differential equations, which we have here analysed in the situations 
where it is possible to decouple a subset of the modes. We have also analysed the large spin 
meson spectra in presence of a magnetic field to realize the well-known Zeeman effect. We 
also found that an external magnetic field enhances the stability of mesons and inhibits 
the dissociation. In appendix D we have presented a model calculation analyzing spinning 
strings in Rindler space to capture key features of the large spin meson dissociation. 

There are several directions to pursue in future work. It is of interest to obtain the low 
spin meson spectra in the low temperature phase {i.e., when chiral symmetry is broken) 
and their dispersion relations. In the present work, we have focussed on the effect of a 
magnetic field. The analysis of the meson spectra and quasinormal modes in the presence 
of an electric field could be pursued in a similar spirit. The effect on the spectra of the 
presence of a chemical potential would also be of great interest. These external fields can 
all be realized as the non-normalizable modes of anti-symmetric fields on the world volume 
of the probe brane, and so we would expect that in all these cases in the high temperature 
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chirally symmetric phase, the scalar and vector fluctuations always remain decoupled, as we 
have seen here. 
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Appendix A 

Here we provide the details of the computation for quadratic action for fluctuation of the 
probe brane. To keep the story general we consider the scalar and vector fluctuation simulta- 
neously since in presence of an external magnetic field they are likely to couple. Parametriz- 
ing the world volume of the DS/DS-hrane by the function To(w), we consider the following 
fluctuation 

t{u) = To{u) + (2W)x(r) , (A-1) 

where To{u) is the classical embedding of the probe brane and {^"-j refer to its worldvolume 
coordinates. 

To proceed let us recall that the DBI lagrangian is given by 



:dbi = e-'^^-det{F[G^, + B^,] + {27ra'Fab)) 

= e-^-det [e^^ + {2^a')E^^ + {2^arE^) 

= e-^^-det (^1 + ^TrM - h^rM^ + ^{TrMf^ + O {M") , (A-2) 



where 



Ml = (2W) {Ej^yE^ + (2W)^ , 

= Grr {daX) {d,X) ■ (A-3) 

Here G^i, and {X'^} refer to the bulk spacetime. For organizational convenience let us also 
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note that 

IIQ4II =diag|(Vi?)-^/'^/^(Vi^)"'^'^i^(V^)~^^'^i^(V^)~'^'^i'} , 

^-det(£;S^) = + ^2) (det||n4||) , (A-4) 

where 5 and A are the symmetric and anti-symmetric part respectively and ||^^4|| denotes 
the diagonal metric for the 5*4. 

For the low temperature background we have 
Grr = {u/Rf f{u) , G,, = Gu = {u/Rf , G^u = (1)"'^' , f{u) = 1 - (^)' 
and for the high temperature background we have 

Gu = {u/Rf f{u) , G,, = Grr = iu/Rf/' , G„„ = (|)"'^' , f{u) = 1 - . 

In what follows we will set i? = 1. 

Now we can determine the effective lagrangian corresponding to small fluctuations. The 
contribution at order {2TTa'y is summarised below 

= ^Grr {S-^GrrT'oiuf - l) e't' ^ -det [e^^) S^' (9„x) {d,x) , 



e-X{u)GrrS''^^-det[E^^)A 



xFab . (A-5) 



The term C^p is the interaction term that couples the scalar and the vector modes. Clearly 
if Tq(m) = then the coupling vanishes. Therefore for the high temperature phase (i.e., 
when Tq{u) — const.) the scalar and vector meson modes are always decoupled. 

Now we determine the contribution coming from the Wess-Zumino term. To this end let 
us first recall that the background has a C5 potential and a C3 potential dual to C5. The 
C5 potential is given by 

C5 = ^ dx° A dx^ A dx"^ A dx^ A dr . (A-6) 

27 



On the worldvolume of the probe brane this induces a 5-form given by 

-^J irQ{u) + 27r a' x)dx° A dx^ A dx'^ A dx^ A du . (A-7) 

Therefore we can immediately observe that for the high temperature phase this cannot 
contribute at the leading order. 

Let us now recall that the Wess-Zumino term is given by 

Swz = yU8 / V A exp (27ra'F + B) . (A-8) 

It is straightforward to see that P[C(5)] A 5(2) vanishes identically. Therefore the total 
contribution coming from the Wess-Zumino term at order {27ra'y is given by 

Swz = ^f^sj (P[C5]AF2AF2-^F^AAAF2AB2^ , with 
F4 = dCs F2 = dA, (A-9) 

We have performed an integration by parts to obtain the second term above. In what follows 
we will focus on the special case of Tq(m) = in which case the first term in eqn. (lA-QO does 
not contribute at the leading order. 

The equations of motion for the gauge fields then decouple from the scalar fluctuation and 
is given by (using R = 1) 



da 



e-^-det (i^l?)^'^"'^^'''^,,,, 



+ 4g;^Be^^''''^^db,Aa' = . (A-10) 



where now a,a',b,b' G {M^''^} |J{m}, on the probe brane worldvolume. It is now possible 
to consistently set all the Aa = where a G S*^. This choice imposes a constraint of the 
following form 

S'^^'d^d.A, = . (A-11) 

This constraint can be trivially satisfied by looking at the gauge field fluctuations indepen- 
dent of the spherical directions (in other words focussing only on the 50(5) singlet states). 

The most general ansatz for the gauge field fluctuations consistent with the previous 
constraint is given by 

_ j^^^\^^-iujt+ikx'^+ik2x'^+ik3X^ 

Au = A^(^)e-*^*+*'=^'+*'=2^'+^'=3x3 _ (A-12) 
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To fix the residual gauge, we further impose Au = 0. This gauge choice gives the following 
constraint 



e-y-det [e^^)s-- [S'\ujA'^ - S'\kA[) - S'\k,A,) - S'\k,A,)] 

+Ag-^B{ujAi + kAt) = . (A-13) 

With this ansatz the equations of motion for the gauge fields are given by 



+ 4g:'BA[ 



-e-M-det [kS^^uAi + kAt) + k2S^\ujA2 + k2At) + k^S^'^uA^ + k^At)] = , 



du 



-e-M -det [k2S^\k2Ai - kA2) + k^S^^k^A^ - kA^) + cuS'^uA^ + kAt)] = , 



du 



e-^/-det 



-e-<^^-det (^1^)522 [kS'\kA2 - k2Ai) + ksS^^k3A2 - k2As) + ujS'\ujA2 + ^2^)] = 
du 



e-^/-det (e2M5335"M'3 



-e-M-det (^1?)533 [kS'\kAs - k,A,) + k2S^^ik2As - k,A2) + ^^"(^^3 + k^At)] = 



(A-14) 



In general, the gauge field fluctuations can have momentum along any of the spatial di- 
rections. Due to the presence of the magnetic field there are two particularly interesting 
cases to consider. The fluctuation field can either oscillate in the direction parallel to the 
magnetic field, or in the direction perpendicular to the magnetic field. We discuss these 
cases below. 

Case 1: Oscillation parallel to the magnetic field 



Let us consider first restricting ourselves to k2 = = k^. Equations in (1A-I4p then reduce 
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to 



du 
du 
du 



-det (^^2^)5**5"" a; + 4.g-^BA^ - e'^W-det {E^f^^S'^kS^^^ujA^ + kAt) = , 



—det 


y ab J 


—det 




—det 





e-M-det - 4:g;^BA[ - e~^J~det (^E^J^S^^ujS'\tuAi + kAt) = , 



- e-M-det [e^^)s'^ {k'S'' + uj'S'') A, = 



-det [E^AS'^^S'^'^A'^ 



along with the following constraint 



e"M-det {k'S'' + iu^S'') A3 = , (A-15) 



e-^^^-det (^l?)^"" [S'\uA[) - S^\kA\)] + Ag;^B{LuA, + kAt) = . (A-16) 

In this case, the longitudinal modes {At,Ai components oscillating along {t,x^} plane) 
and the transverse modes (^42,^43 components oscillating along {t,x^} plane) clearly decou- 
ple. It is useful to note that A2, A^ transform as vectors and {At, Ai} transforms as a scalar 
under the SO (2) rotation group in the plane perpendicular to the magnetic field. 

Case 2: Oscillation perpendicular to the magnetic field 



Now we restrict ourselves to /c = = /c2. In this case the equations of motion in (IA-14P 
become 



du 
du 
du 
du 



e-M-det (^E^^y^'S^^A't + Ag;^BA[ - e-M-det (^E^^y'^ksS^^uAs + k^At) = 



e-^^-det - Ag^'BA't - e^^-^-det (^l?)^" {klS'^ + u'S'') A, = 

■det lE'rAs^^S'''' A' 







— e 



-det (^^3 + ksAt) = , 



-det(E2M5225"M'2 



-det (klS^^ + cu^S'') = , (A-17) 



along with the constraint 



e~y -det [S'\ujA't) - S^^ik^A',)] + Ag;'BiLuA^ + kAt) = 

30 



In this case only the y42-transverse mode decouple from the rest of the gauge field fluctu- 
ations. The longitudinal modes {At,As components oscillating along {t,x^} plane) and the 
transverse mode {Ai component oscillating along {t, x^} plane) are in general coupled. 



Appendix B 

Here we analyze the gauge field fluctuations when there is no magnetic field. We provide 
an analytical derivation for the dispersion relation of the lowest hydrodynamic mode of the 
longitudinal oscillation, which has been obtained by other methods in refs. [5U] . 

The longitudinal mode is defined to be S{u) = uAi + kAf. Using this definition, the 
constraint equation in (1A-16P and the equations of motion in (lA-lSp (after setting B = 0) 
we obtain the equation of motion for the longitudinal mode to be 

£iu){u'-ef{u))R' £'iu){5f{u)u' + 2uf'{u)u^-5k^f{ur) 

u^f{uf 2uf{u){u^-k^f{u)) +^ W u, 

-) 

With the following variable change 



/W = l-(^y. (B-l) 



x = —, (B-2) 

Ut 

the equation of motion for the longitudinal mode becomes 

(PS_ 5 (1 - x-^) {uj^ - P (1 - x-3)) + 6u;^x-^ dS R^oo^ - k^l - x-^)) 

^ 2x(l-x-3) (cu2 - p (1 - x-3)) + ~ ■ 

(B-3) 

The near-horizon limit u — > is now achieved by taking the x 1 limit. In this limit, 
the equation of motion takes the form 

d^S _^ ^ dS ^ R^u^ ^^^^ _ ^ 
dx"^ X — 1 dx 9ut{x — ly 

The general solution of this equation is given by a linear combination of the incoming and 
the outgoing modes 



£(x) = Cicos(cI;log(x- 1)) +C2sin(tZ'log(x- 1)) , Cu = J — = — 77; • (B-5) 

y 9ut 47rT 

The incoming boundary condition singles out the solution with Ci = — C2. Near the horizon 
S{x) is therefore obtained to be 

£{x) = CiExp {~iCd log{x - 1)) , (B-6) 
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where is an yet undetermined constant. 

Let us also define k = k/{A7TT). Now, for small enough uj and k (such that ^ 1 and 
k <^ 1), we can ignore the last term in equation (lB-3p . It turns out that in this limit the 
equation of motion for the longitudinal mode is also exactly solvable. 



S{x) = M2 + -Ml 



X 



3/2 



+ uj"^ log 



X' 



3/^ + 1 

3/2 _ 1 



(B-7) 



where Mi and M2 are constants of integration. Near the boundary, i.e., x ^ 00, this reduces 
to 



S{x) = M2 - -Ml (P -u^^ ^ 
3 V 



3/2 



X 



(B-8) 



Therefore normalizability of S{x) forces us to impose M2 = 0. The quasi- normal modes are 
therefore obtained to be the solution of this constraint. 



On the other hand, in the vicinity of the horizon the solution in (IB-7P takes the following 
form 

'4^ 



Mia;^log(a; - 1) + 



S{x) = - ( -2MiP + 3M2 + Miu^ log 

3 Y \ 3 

^Mi ( 6P + u^) (x - 1) - -i-Mi ( I8OP + u^){x- 1? 
6 V / 144 



(B-9) 



Now, for sufficiently small values of a;, we can expand the solution obtained in (1B-6P to 
get 

1 



Six) =Cl- tCLO\og{x - 1) - ^ClCu' i\og{x - l)y 



(B-10) 



Comparing equation flB-9p and fIB-lOp we get 

1 



-2Mik^ + 3M2 + Mit:;^ 
3 Y \ 3 



Setting M2 = 0, the solution of this equation is given by 

Mi = ^, cD = -2^P + ... 
uj 



Cr 



(B-11) 



(B-12) 



Restating this result in dimensionful parameters, we get the lowest hydrodynamic quasinor- 
mal frequency 

1 



u = —iDnk , D} 



2'kT 



(B-13) 
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where Dji is the R-charge diffusion constant. This value for the diffusion constant was 
obtained numerically by studying the spectral functions of holographic flavours in ref. [15] 
and was also confirmed in ref. [30] by numerically obtaining the lowest quasinormal mode 
of the longitudinal oscillation. 

In the vanishing magnetic field limit, the transverse modes obey the following equations 
(obtained from eqn. (lA-lSp ) 

A,{u)iu;'-k'fiu))R' , A',{u){5f{u) + 2uf'{u)) 



+ Al{u) = 0. (B-14) 



u^fiuy 2uf{u) ' 

A,{u){uj'-kyXu))R' , Ai,iu){5f{u) + 2uf'{u)) 



m3/(m)2 2uf{u) 

Clearly the A2 and A^ meson spectra are degenerate. For convenience we again use the 
variable x = u/ut- The solution close to the horizon is similar to the one obtained for the 
longitudinal mode (given in equation ( lB-61) ) and is given by 

A2{x) = CtExp {~iCd log(x - 1)) , (B-15) 

where Ct is an yet undetermined constant. At small values of u and k (compared to the 

background temperature T) the solution for the transverse modes can be obtained to be 

-3,3/2 ^ ^- 



A2{x) =N2-^Ni\o^ 



X' 



3/2 _ 1 



+ . . . , (B-16) 



where A^i and N2 are constants of integration. 

Now expanding the solution in eqn. flB-161) near the horizon and matching with the 
solution in eqn. flB-lSp gives 

N2-^N,\og(^^^ =Ct, 

-Ni = -iCru . (B-17) 
3 

Normalizability requires us to set N2 = 0. However, we do not get a consistent solution for 
A'^i from the conditions above. From the equations of motion in ( ]B-14[) . it can be checked 
numerically (see also ref. [3t)J) that the normalizability condition cannot be satisfied in the 
hydrodynamic limit (namely when uj <^ 1 and k <^ 1). Therefore we can conclude that the 
transverse mode does not have any solution compatible with the limit -C 1 and A; <^ 1 in 
accordance with the results obtained in, e.g. ref. [55] . 



Appendix C 

Here we outline the general variable changes to obtain the corresponding Schrodinger equa- 
tion from the original equation for fluctuation modes. Without any loss of generality we 
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can write the equation for any fluctuation mode as follows 



ai(w)/v(w) + a2iu)f'{u) + uj'^a3{u)f^{u) 







(C-1) 



where, ai{u), 02 (w) and a^lu) are known functions and prime denotes derivative with respect 
to u, fx{u) denotes the radial profile of the fluctuation mode and uj is the corresponding 
oscillatory frequency (in the finite temperature case this is the quasinormal mode). Let us 
rewrite f^iu) = a{u)g{u) with 



a\u) 



CTiU 



ai{u) 2 



ai[u) 
03 (m) 



du 



ai{u) 



(C-2) 



The Schrodinger equation is then obtained to be 




{dug{uj) + uj'^g{u) - Vs{u)g{u) = , where 



A{u) + 



B u) 



B{u) 



+ 2'^"l + Bin) 



1/2 



ai[u) \ai[u) J 

A more conventional form of the Schrodinger equation can be obtained after changing vari- 
ables to "tortoise" coordinate du = B{u)du] the horizon is then located at n ^ 00. The 
main goal of this exercise is to obtain the potential Vs{u) and extract qualitative features 
of the meson spectrum. 

Here we have explicitly assumed that the fluctuation does not have any momentum mode. 
The fluctuation equation for the meson having a momentum can also be recast in the form 
of a Schrodinger equation. The changes of variables are exactly similar as already mentioned 
in eqn. (IC-Sp . but the effective potential receives a positive contribution coming from the 
momentum 



as 



(C-4) 



where the left hand side of the equation ( IC-ip is now accompanied by another term of the 
form —\a4{u)\k'^fx, the relative negative sign is due to the fact that /c-corresponds to a 
spatial oscillation. 



Appendix D 

Here we would like to solve a toy problem to study the effect of the magnetic field on meson 
dissociation. We implicitly assume the framework where we can have two equal charges at 



34 



the two ends of the string and therefore analyze the symmetric configuration only. Rindler 
spacetime provides an useful arena where many qualitative (if not quantitative) features 
of an event horizon can be realized in a simple set up. Such study has been previously 
carried out in ref. [28] without any background field. The main result of this exercise is that 
the high spin meson dissociates once a critical value of the acceleration is reached. This 
acceleration is in turn determined by the angular momentum of the meson. 

There are two ways to physically interpret the results of such pursuits (e.g., see ref. [28] and 
ref. [27]). We consider accelerating the rotating string in one of the space-like directions and 
therefore a Rindler horizon forms in the four dimensional spacetime. Alternatively, we can 
consider accelerating the string in the holographic direction and therefore a bulk spacetime 
horizon forms which sets the temperature of the dual gauge theory. To begin with we take 
the former point of view. 

The metric for the Rindler space can be written as 

ds^ = -('K'dTf + d^" + dp" + p2d02 ^ 

where {p, 0}-plane represents the plane where we would rotate the string. To add a back- 
ground magnetic field we consider similar potential given before. Now the ansatz for the 
string in terms of its worldvolume coordinates {r, cr} is given by 

7] = r , i = ^(o-) , p = p{a) , (j) = ujT . (D-2) 

This ansatz implies that the direction of acceleration of the string is orthogonal to the direc- 
tion of its rotation. We can readily see that there are two conserved charges associated with 
the string, namely the angular momentum J, which is associated with the rotation along 
0-direction and the boost charge, which is associated with the translation in ?7-direction. 

The Nambu-Goto action for the string now will be accompanied by a boundary term 
exactly similar to the term ASb given in eqn. (fT6|) . Therefore the full action is given by 

S = ^ [ rfrdaV(e'«'-^V)(f' + P") + A5b . (D-3) 
Zna J 

Varying this action we can obtain the equation of motion for the string profile and the 
boundary conditions we should impose. In absence of any external field this exercise has 
been carried out in details in ref. [28], which we do not repeat here. The presence of the 
magnetic field enforces us to impose the following boundary condition 



- ^ 



= -Bpuo , (D-4) 

±£/2 
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where £ is the distance between the two end points of the string. The other boundary 
term corresponding to the variation of C, is satisfied imposing the Dirichlet boundary con- 
dition 6^ = 0. This corresponds to a situation where the end points move with a constant 
acceleration a = ^{±i/2)~^. 

Now choosing a gauge p = a we can 

solvfl 

the equation of motion to obtain the profile 

for the string given by 



As the notation suggests the parameter J represents the angular momentum of the meson. 

As we have seen before, the presence of the background magnetic field does not change 
the equations of motion for the Nambu-Goto string. Therefore the functional form of the 
string profile given in equation (]D-5I) is entirely determined by the Rindler geometry [28] . We 
notice that the presence of the magnetic field does induce an effective length £eff = i'^{l + B^) 
by changing the boundary condition of the string end-points. 

A non-zero value of B brings about two key changes in the evaluation of the angular 
momentum. The profile of the string modifies in order to satisfy the boundary condition 
in eqn. (1D-4I) and therefore the effective length ies arises in the formula in ( ]D-5I) . There is 
also an additional explicit contribution AJb of the form shown in eqn. f|T8l) . which finally 
gives a functional factor of tan^^(l/i?) in the expression for the angular momentum. 

It can be shown that the explicit contribution of the of the magnetic field to the angular 
momentum, denoted by AJb, is actually cancelled by a term coming from integrating over 
the string profile. The net resulting angular momentum therefore can be lower than the 
angular momentum of the string in absence of the field. This is pictorially summarised in 
figure [151 

To find the maximum acceleration we identify a^^ = and look for minima of the 

right hand side of the equation as a function of the parameter C. For now we content 
ourselves with only positive values of B, i.e. when the string angular momentum and the 
background field are parallel. It is a straightforward exercise to show that such minima 
corresponds to values of C which satisfy the following equation (obtained from taking the 

-•^^The analytical solution of the equations of motion obtained from the Nambu-Goto action remains the 



same as in ref. [28]. We refer to ref. [28] for further details about solving the equations of motion. We use 
the same analytical solution and impose the boundary condition given in (|D-4p to obtain the results in eqn. 




(D-5) 



jEll). 
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0.44 




Figure 15: The total angular momentum of the string in units of (27ra') in the background field. 
The horizontal dashed line represents the value of angular momentum of the string when there is 
no background magnetic field. We observe that the B-Held initially reduces and then increases the 
angular momentum creating a local minima. However, we will observe that the stability of the 
meson is always enhanced by the 5-field. 



first derivative of and setting it equal to zero at p = ±^/2) 
a; = mtanh^ — ^ , where 

CI ( \ 

X = . , m = arcsm 



K V2WJ \y/l + B 



arctan | — 
B 



-1/2 



(D-6) 



Clearly the parameter B generates a set of such values of C which satisfy the above equation 
and therefore promotes the maximum value of acceleration to a function of the magnetic 
field. Differentiating the eqn. in the first line of flD-6|) and using identities for the trigono- 
metric hyperbolic functions we obtain that the family of roots of the equation can be simply 
given by the relation x = am where a = 0.834 is a constant. 

Using these roots the maximum acceleration can be obtained as a function of the back- 
ground magnetic field which is given by 



oc ^/arctan I -—j j ^ . (D-7) 

V / arcsin [^^jig^j 

The functional dependence is pictorially shown in figure [161 We observe that the maximum 
acceleration grows almost linearly with the background field. 

Alternatively we could imagine the Rindler space to be the near-horizon approximation of 
a blackhole background and the acceleration is along the holographic coordinate (we make 
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the identification that ^ = m, where u is the radial coordinate in the holographic set up). 
The parameter k (in eqn. fID-ll) ) in this case represents the surface gravity and sets the 
temperature of the dual gauge theory to be T = k/{2t[). 

Now we imagine the string end points to sit on the flavour brane at some value of = Uq, 
which fixes the constituent quark mass; and we consider a similar ansatz for the string profile 
as presented in eqn. (lD-2p . The interpretation of the set up is different from the earlier one; 
however this brings about no change as far as the mathematics is concerned. We therefore 
again obtain the string profile as given by eqn. (]D-5|) . 

However in this case the physical parameter that captures the dissociation of the meson 
is the size I (which sets the angular momentum J via the relation in eqn. (ID-SP ) of it. The 
length is obtained to be 



2C arccosh [kUq/C) 



arcsin (l/Vl + B^^ ^ . (D-8) 



We would hope to see that for a given value of Uq, there exists a maximum admissible size 
of the spinning meson beyond which it dissociates. 

We can treat as before the parameter C to be independent with respect to which we will 
consider maximizing the length function i. From eqn. ( ID-Sj) it is clear that the external 
magnetic field does not play any role to determine the value of C corresponding to the 
maximum of i; however it determines the the function £max(-B) and therefore also Jmax(-B) 
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to be given by 
1 



oc 



arcsm (i/Vi + 52 



, Jmax OC arctan ( — 



arcsm (1/Vi + 52 



:d-9) 



where the constants of proportionahty depends on the background temperature T and the 
parameter Uq- 




Figure 17: We plot the two functions given in eqn. ()D-9p . The constants of proportionality have 
not been taken into account here. 



For illustrative purpose we have plotted the functional behaviour of £max and Jmax in figure 



17(a) and 17(b) respectively. We see that the maximum length of the meson saturates an 
upper bound at high enough values of the magnetic field; however no such saturation is 
present in the maximum angular momentum of the meson. 

For the background considered in the main text the Rindler spacetime does emerge when 
we zoom in the region close to the horizon [T3] . The toy model here can therefore be 
identified with a corresponding study of high spin mesons in the overheated phase. So if we 
want to make a direct connection between the Rindler background dynamics to the specific 
system studied in the text, we should remember that such a background appears precisely 
in the over-heated phase of the system, where any small fiuctuation is likely to destroy the 
meson and drive the over-heated phase to the melted phase. 
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